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Abstract. We study the geodesic X-ray transform Ir of tensor fields on a compact Riemannian 
manifold M with non-necessarily convex boundary and with possible conjugate points. We assume 
that Ir is known for geodesies belonging to an open set T with endpoints on the boundary. We 
prove generic s-injectivity and a stability estimate under some topological assumptions and under 
the condition that for any (x,^) £ T*M, there is a geodesic in F through x normal to ^ without 
conjugate points. 



Let (M, dM) be a smooth compact manifold with boundary, and let g G C^{M) be a Riemannian 
metric on it. We can always assume that (M, dM) is equipped with a real analytic atlas, while 
dM and g may or may not be analytic. We define the geodesic X-ray transform / of symmetric 
2-tensor fields by 



where [0, l-y] 9 1 1-^ 7 is any geodesic with endpoints on dM parameterized by its arc-length. Above, 
ifjO"^) is the action of / on the vector 6, that in local coordinates is given by fijO^OK The purpose 
of this work is to study the injectivity, up to potential fields, and stability estimates for / restricted 
to certain subsets T (that we call Iy)-, and for manifolds with possible conjugate points. We require 
however that the geodesies in V do not have conjugate points. We also require that T is an open sets 
of geodesies such that the collection of their conormal bundles covers T*M. This guarantees that 
Ir resolves the singularities. The main results are injectivity up to a potential field and stability 
for generic metrics, and in particular for real analytic ones. 

We are motivated here by the boundary rigidity problem: to recover g, up to an isometry 
leaving dM fixed, from knowledge of the boundary distance function p{x, y) for a subset of pairs 
€ dM X dM, see e.g., [Mil IShl( ICDSt [SU41 IPU| . In presence of conjugate points, one should 
study instead the lens rigidity problem: a recovery of g from its scattering relation restricted to a 
subset. Then I-p is the linearization of those problems for an appropriate F. Since we want to trace 
the dependence of /r on perturbations of the metric, it is more convenient to work with open F's that 
have dimension larger than n, if n > 3, making the linear inverse problem formally overdetermined. 
One can use the same method to study restrictions of I on n dimensional subvarieties but this is 
behind the scope of this work. 

Any symmetric 2-tensor field / can be written as an orthogonal sum of a solenoidal part and 
a potential one dv, where f = on dM, and d stands for the symmetric differential of the 1-form 
V, see Section [21 Then I{dv){'y) = for any geodesic 7 with endpoints on dM. We say that Ip is 
s-injective, if /p/ = implies f = dv with v = on dM, or, equivalently, / = This problem 
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has been studied before for simple manifolds with boundary, i.e., under the assumption that dM 
is strictly convex, and there are no conjugate points in M (then M is diffeomorphic to a ball). 
The book |Shl] contains the main results up to 1994 on the integral geometry problem considered 
in this paper. Some recent results include |Sh2| . |Ch] . [SU3| . [D] . |Pe| . |SSU| . |ShU] . In the two 
dimensional case, following the method used in |PUj to solve the boundary rigidity problem for 
simple 2D manifolds, injectivity of the solenoidal part of the tensor field of order two was proven in 
|Sh3| . In |SU4j , we considered / on all geodesies and proved that the set of simple metrics on a fixed 
manifold for which / is s-injective is generic in C^{M), k ^ 1. Previous results include s-injectivity 
for simple manifolds with curvature satisfying some explicit upper bounds |ShH ISh21 IPe| . A recent 
result by Dairbekov [DJ proves s-injectivity for non-trapping manifolds (not-necessarily convex) 
satisfying similar bounds, that in particular prevent the existence of conjugate points. 

Fix another compact manifold Mi with boundary such that M™* D M, where M™* stands for 
the interior of Mi. Such a manifold is easy to construct in local charts, then glued together. 

Definition 1. We say that the C^{M) (or analytic) metric g on M is regular, if g has a (or 

analytic, respectively) extension on Mi, such that for any (x,^) E T*M there exists 9 G T^M \ 
with 0) = such that there is a geodesic segment ^^^q through (x, 6) such that 

(a) the endpoints of jx,e o-fc in M^^ \ M . 

(h) there are no conjugate points on ^xfi- 
Any geodesic satisfying (a), (b) is called a simple geodesic. 

Note that we allow the geodesies in V to self-intersect. 

Since we do not assume that M is convex, given {x,6) there might be two or more geodesic 
segments 7j issued from (x, 0) such that 7^ n M have different numbers of connected components. 
Some of them might be simple, others might be not. For example for a kidney-shaped domain and 
a fixed (x, 9) we may have such segments so that the intersection with M has only one, or two 
connected components. Depending on which point in T*M we target to recover the singularities, 
we may need the first, or the second extension. So simple geodesic segments through some x (that 
we call simple geodesies through x) are uniquely determined by an initial point x and a direction 
9 and its endpoints. In case of simple manifolds, the endpoints (of the only connected component 
in M, unless the geodesies does not intersect M) are not needed, they are a function of {x,9). 
Another way to determine a simple geodesic is by parametrizing it with (x, rj) £ T(Af™* \ M), such 
that exp^ 7] G M'^^ \ M then 

(2) 7,,^ = {expjtry),0<t< 1}. 

This parametrization induces a topology on the set F of simple geodesies through points of M}"*. 

Definition 2. The set F of geodesies is called complete, if 

(a) V(x,^) G T*M there exists a simple geodesic 7 G F through x such that 7 is normal to ^ at 

X. 

(h) T is open. 

In other words, a regular metric g'lsa, metric for which a complete set of geodesies exists. Another 
way to express (a) is to say that 

(3) A^*F := {iV*7; 7 G F} D T*M, 

where A^*7 stands for the conormal bundle of 7. 

We always assume that all tensor fields defined in M are extended as to Afi \ M. Notice 
that If does not change if we replace M by another manifold Mi 12 close enough to M such that 
M C Mii2 C Ml but keep / supported in M. Therefore, assuming that M has an analytic 
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structure as before, we can always extend M a bit to make the boundary analytic and this would 
keep (M, dM, g) regular. Then s-injectivity in the extended M would imply the same in the original 
M, see |SU41 Prop. 4.3]. So from now on, we will assume that (M, dM) is analytic but g does not 
need to be analytic. To define correctly a norm in C^{M), respectively C^{Mi), we fix a finite 
analytic atlas. 

The motivation behind Definitions [U [2] is the following: if g is regular, and V is any complete set 
of geodesies, we will show that Ip/ = implies that /'^ G C'(M), where / = l{k) — > oo, as A; — > oo, 
in other words, the so restricted X-ray transform resolves the singularities. 

The condition of g being regular is an open one for g S C^{M), i.e., it defines an open set. Any 
simple metric on M is regular but the class of regular metrics is substantially larger if dim M > 3 
and allows manifolds not necessarily diffeomorphic to a ball. For regular metrics on M, we do not 
impose convexity assumptions on the boundary; conjugate points are allowed as far as the metric 
is regular; M does not need to be non-trapping. In two dimensions, a regular metric can not have 
conjugate points in M but the class is still larger than that of simple metrics because we do not 
require strong convexity of dM. 

Example 1. To construct a manifold with a regular metric g that has conjugate points, let us start 
with a manifold of dimension at least three with at least one pair of conjugate points u and t; on a 
geodesic [a, 6] 3 lif)- We assume that 7 is non-selfintersecting. Then we will construct M as 
a tubular neighborhood of 7. For any xq S 7, define = exp^^^jf; (f,7(xo)) = 0, |f| < e}, and 
M := VJxQ^^SxQ with e <C 1. Then there are no conjugate points along the geodesies that can be 
loosely described as those "almost perpendicular" to 7 but not necessarily intersecting 7; and the 
union of their conormal bundles covers T*M. More precisely, fix x S M, then x G Sxq for some 
xq G 7. Let 7^ ^ E T*M. Then there exists 7^ v S TxM that is both tangent to S^o s-^id normal 
to ^. The geodesic through [x, v) is then a simple one for e ^ 1, and the latter can be chosen in a 
uniform way independent of x. To obtain a smooth boundary, one can perturb M so that the new 
manifold is still regular. 

Example 2. This is similar to the example above but we consider a neighborhood of a periodic 
trajectory. Let M = {^{x^)"^ + (x^)^ < l} x be the interior of the torus in R^, with the fiat 
metric (dx^)^ -|- (dx^)^ + d9'^, where 6 is the natural coordinate on with period 27r. All geodesies 
perpendicular to ^ = const, are periodic. All geodesies perpendicular to them have lengths not 
exceeding 2 and their conormal bundles cover the entire T*M (to cover the boundary points, we 
do need to extend the geodesies in a neighborhood of M). Then M is a regular manifold that is 
trapping, and one can easily show that a small enough perturbation of M is also regular, and may 
still be trapping. 

The examples above are partial cases of a more general one. Let (M' , dM') be a simple compact 
Riemannian manifold with boundary with dimM' > 2, and let M" be a Riemannian compact 
manifold with or without boundary. Let M be a small enough perturbation of M' x M" . Then M 
is regular. 

Let g he a fixed regular metric on M. The property of 7 being simple is stable under small 
perturbations. The parametrization by (x, rj) as in ([2]) clearly has two more dimensions that what 
is needed to determine uniquely 7|m- Indeed, a parallel transport of (x,r/) along 'yx,ri, close enough 
to X, will not change 7|mi similarly, we can replace ??by(l-|-e)7/, 

We assume throughout this paper that M satisfies the following. 

Topological Condition: Any path in M connecting two boundary points is homotopic to a 
polygon ci U 71 U C2 U 72 U • • • U 7fc U c^+i with the properties: 
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(i) Cj are paths on dM; 

(ii) For any j, 7j = 7j|M for some 7j S F; 7^ lie in M™* with the exception of its endpoints and 
is transversal to dM at both ends. 

Theorem 1. Let g be an analytic, regular metric on M . Let T be a complete complex of geodesies. 
Then /r is s-injective. 

The proof is based on using analytic pseudo-differential calculus, see [Sj, ITrej . This has been 
used before in integral geometry, see e.g., [BQ Q], see also |SU4j . 



To formulate a stability estimate, we will parametrize the simple geodesies in a way that will 
remove the extra two dimensions. Let Hm be a finite collection of smooth hypersurfaces in M™*. 
Let Hm be an open subset of {{z,9) G SMi; z G Hm,0 TzHm}, and let ±l^{z,9) > be two 
continuous functions. Let T(TCm) be the set of geodesies 

(4) r(w„) = {j,,e{t); im(.z,e) < t < i+{z,d), {z,e) g n^} , 

that, depending on the context, is considered either as a family of curves, or as a point set. We 
also assume that each 7 G T(7im) is a simple geodesic. 

If g is simple, then one can take a single H = dMi with 1^=0 and an appropriate l~^{z,9). If 
g is regular only, and T is any complete set of geodesies, then any small enough neighborhood of a 
simple geodesic in T has the properties listed above and by a compactness argument on can choose 
a finite complete set of such r('Hm)'s, that is included in the original F, see Lemma [H 

Given 7i = {Tim} as above, we consider an open set TC' = {Ti'^}, such that H'm ^ Hm, and let 
T{7i'm) be the associated set of geodesies defined as in (jH), with the same l^. Set T{7i) = UT{7im)i 

T{n') = urin'm). 

The restriction 7 G T{Ti.'^) C TiTim) can be modeled by introducing a weight function Om in 
Hm, such that am = 1 on Ti'm, and Om = otherwise. More generally, we allow am to be smooth 
but still supported in Tim- We then write a = {am}, and we say that a G C^{H), if am G C^{7im), 
Vm. 

We consider La^ = amL, or more precisely, in the coordinates {z, 9) G TLm, 

(5) L^^f = am{z,9) {f{j,,e),ile)dt, {z,9) ^Hm- 

Jo 

Next, we set 

(6) Ia = {IaJ, N^^=LlJ^^=L*\am?L, N^ = J2n^^, 

where the adjoint is taken w.r.t. the measure d^ := \{i'{z),9)\dSzd9 on Tim, dSzd9 being the 
induced measure on SM, and being a unit normal to Hm.- 

S-injectivity of is equivalent to s-injectivity for /„, which in turn is equivalent to s-injectivity 
of / restricted to suppa, see Lemma [2j The space is defined in Section [2l see ([8]). 

Theorem 2. 

(a) Let g = go ^ , k ^ I be regular, and let 7i' <^ 7i be as above with r{7i') complete. Fix 
Oi = {am} G with Ti'm C suppa^ C Tim- Then if L^ is s-injective, we have 

(7) \\r\\LHM)<C\\No^f\\f,,^,,^y 

(b) Assume that a = ag in (a) depends on g ^ , so that C^{Mi) 3 g ^ C\7i) 3 ag is 
continuous with Z » 1, A; » 1. Assume that Lg^^ag^ is s-injective. Then estimate ^ remains true 
for g in a small enough neighborhood of go in C^{Mi) with a uniform constant C > 0. 
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In particular, Theorem [2] proves a locally uniform stability estimate for the class of non-trapping 
manifolds considered in [P] . 

Theorems disallow us to formulate generic uniqueness results. One of them is formulated below. 
Given a family of metrics Q C C^{Mi), and Ug C T(M™* \ M), depending on the metric g ^ G, 
we say that Ug depends continuously on g, if for any go € G, and any compact K C Ug^^, we have 
K C C/™* for fl' in a small enough neighborhood of go in C^. In the next theorem, we take Ug = T^, 
that is identified with the corresponding set of {x,ri) as in ([2]). 

Theorem 3. Let Q C C^{Mi) he an open set of regular metrics on M , and let for each g ^ Q, Tg 
he a complete set of geodesies related to g and continuously depending on g. Then for k ^ 0, there 
is an open and dense subset Qs of Q, such that the corresponding X-ray transform Iy^ is s-injective. 

Of course, the set Gs includes all real analytic metrics in Q. 

Corollary 1. Let TZ{M) be the set of all regular metrics on M equipped with the C^(Mi) 
topology. Then for k ^ \, the subset of metrics for which the X-ray transform I over all simple 
geodesies is s-injective, is open and dense in TZ{M). 

The results above extend the generic results in [SU4] . see also |SU3] . in several directions: the 
topology of M may not be trivial, we allow conjugate points but we use only geodesies without 
conjugate points; the boundary does not need to be convex; and we use incomplete data, i.e., we 
use integrals over subsets of geodesies only. 

In Section [H we discuss versions of those results for the X-ray transform of vector fields and 
functions, where the proofs can be simplified. Our results remain true for tensors of any order m, 
the necessary modifications are addressed in the key points of our exposition. To keep the paper 
readable, we restrict ourselves to orders m = 2, 1, 0. 

2. Preliminaries 

We say that / is analytic in some subset U of an analytic manifold, not necessarily open, if / can 
be extended analytically to some open set containing U . Then we write / G A{U). Let g G C^{M), 
k ^ 2 or g £ A{M) be a Riemannian metric in M. We work with symmetric 2-tensors / = {fij} 
and with 1-tenors / differential forms vj (the notation here and below is in any local coordinates) . We 
use freely the Einstein summation convention and the convention for raising and lowering indices. 
We think of fij and f^^ = fkig^^g^-' as different representations of the same tensor. If ^ is a covector 
at X, then its components are denoted by while is defined as = g^-'^j- Next, we denote 
I^P = ^i^*, similarly for vectors that we usually denote by 9. If 9i, 62 are two vectors, then {9i, 62) 
is their inner product. If ^ is a covector, and is a vector, then 9) stands for £,(9). This notation 
choice is partly justified by identifying ^ with a vector, as above. 

The geodesies of g can be also viewed as the x-projections of the bicharacteristics of the Hamil- 
tonian Eg{x,£) = \g^^ {x)£,i£j. The energy level Eg = 1/2 corresponds to parametrization with 
the arc-length parameter. For any geodesic 7, we have f^^{x)£i£^j = fij{'~f{x))'j^{t)'y^ (t), where 
(x,^) = is the bicharacteristic with x-projection equal to 7. 

2.1. Semigeodesic coordinates near a simple geodesic and boundary normal coordi- 
nates. Let [l~,l~^] B t 7a;o,6»o(*) be a simple geodesic through xq = 7x0,6*0 (0) S Mi with 
^'o £ SxqMi. The map t9 1— > exp^g{t9) is a local diffeomorphism for 9 close enough to 6*0 and 
t G [^1^^] by our simplicity assumption but may not be a global one, since 7a;o,6»o self- 
intersect. On the other hand, there can be finitely many intersections only and we can assume 
that each subsequent intersection happens on a different copy of M. In other words, we think of 
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7o as belonging to a new manifold that is a small enough neighborhood of 70, and there are no 
self-intersections there. The local charts of that manifold are defined through the exponential map 
above. Therefore, when working near jxo,9o '^^ assume that 'Jxofio does not intersect itself. We 
will use this in the proof of Proposition [2j Then one can choose a neighborhood U of 70 and normal 
coordinates centered at xq there, denoted by x again, such that the radial lines t t9, 6 = const., 
are geodesies. If 5 £ C^, then we lose two derivatives and the new metric is in C^~'^; if g is analytic 
near 70, then the coordinate change can be chosen to be analytic, as well. 

If in the situation above, let xq ^ M, and moreover, assume that the part of 7x0,00 corresponding 
to t < is still outside M. Then, one can consider {9,t) as polar coordinates on Tx^M. Considering 
them as Cartesian coordinates there, see also |SU31 sec. 9], one gets coordinates near jxoM 

so that the latter is given by {(0, . . . ,0, i), < t < l^}, gin = Sin, and = r"„ = 0, Vz. Given 
X G R", we write x' = (x^, . . . , x"~^). Moreover, the lines x' = const., \x'\ <C 1, x" = t G [0, l^] are 
geodesies in T, as well. We will call those coordinates semigeodesic coordinates near 7xo,6»o- 

We will often use boundary normal (semi-geodesic) coordinates {x',x^) near a boundary point. 
If x' G R"^^ are local coordinates on dM, and i^(x') is the interior unit normal, for p £ M close 
enough to dM, they are defined by exp^^/ x"z/ = p. Then x" = defines dM, x" > in M, 
x" = dist(x, dM). The metric g in those coordinates again satisfies gin = Sin, and FJ^^ = F"^ = 0, 
Vi. We also use the convention that all Greek indices take values from 1 to n — 1. In fact, the 
semigeodesic coordinates in the previous paragraph are boundary normal coordinates to a small 
part of the geodesic ball centered at xq = 7a;o,0o(O) with radius e, < e ^ 1. 

2.2. Integral representation of the normal operator. We define the space of symmetric 
tensors / with inner product 

if,h)= [ {f,h)idetg)'/^dx, 

where, in local coordinates, {f,h) = fijff^ . Similarly, we define the space of 1-tensors (vector 
fields, that we identify with 1-forms) and the LP' space of functions in M. Also, we will work in 
Sobolev spaces of 2-tensors, 1-forms and functions. In order to keep the notation simple, we 
will use the same notation (or H'^) for all those spaces and it will be clear from the context 
which one we mean. 

In the fixed finite atlas on M, extended to Mi, the norms ||/||(7fe and the norms below are 
correctly defined. In the proof, we will work in finitely many coordinate charts because of the 
compactness of M, and this justifies the equivalence of the correspondent and H'^ norms. 

We define the Hilbert space H'^{Mi) used in Theorem [2] as in [STJ3l ISTli] . Let x = (x',x") be 
local coordinates in a neighborhood U of a point on dM such that x" = defines dM. Then we 
set 




This can be extended to a small enough neighborhood V of dM contained in Mi. Then we set 

n 

(8) ll/lli/2(Mi) = X] \\^xjf\\m(V) + ll/lli?i(Mi)- 

The space if^(Mi) has the property that for each / G H^{M) (extended as zero outside M), we 
have Nf G H^{Mi). This is not true if we replace H'^{Mi) by H'^{Mi). 

Lemma 1. Let Tg and Q he as in Theorem Then for k ^ 1, for any go £ G, there exist 
TL' = {H'ni} = {Hm} such that T{7i) (s F^p, and TL' , TL satisfy the assumptions of Theorem\^ 
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Moreover, 7i' and Ti satisfy the assumptions of Theorem{Efor g in a small enough neighborhood of 
go in C^. 

Proof. Fix go (z G first. Given {xq,(^q) G T*M, there is a simple geodesic 7 : Mi in 

Tgg through xq normal to ^0 at xq. Choose a small enough hypersurface H through xq transversal 
to 7 G Tgg, and local coordinates near xq as in Section [2.11 above, so that xq = 0, H is given by 
= 0, 7(0) = (0, . . . , 0, 1). Then one can set Hq = {x; x" = 0; \x'\ < e} x {9; \0'\ < e}, and Ti'^ 
is defined in the same way by replacing e by e/2. We define T{7{o) as in ([4]) with l^{z,0) = . 
Then the properties required for TIq, including the simplicity assumption are satisfied when < 
e <C 1. Choose such an e, and replace it with a smaller one so that those properties are preserved 
under a small perturbation of g. Any point in SM close enough to (xo,^o) still has a geodesic in 
r('HQ) normal to it. By a compactness argument, one can find a finite number of so that the 
corresponding T{'H') = ur('H^) is complete. 

The continuity property of Tg w.r.t. g guarantees that the construction above is stable under a 
small perturbation of g. □ 

Similarly to |SU3| . one can see that the map la^ : L'^{M) L'^{7im, d/x) defined by ([5]) is 
bounded, and therefore the normal operator defined in ([6|) is a well defined bounded operator 
on L'^(M). Applying the same argument to Mi, we see that : M Mi is also bounded. By 
|SU3j . at least when / is supported in the local chart near xq = above, and x is close enough to 



Jo Js^n 

where lamCa^j^*)!^ = \oimix,0)\'^ + \amix, —9)\'^ , and am is the extension of am as constant along 
the geodesic through (x, 0) E Tim', and equal to for all other points not covered by such geodesies. 
Formula ([9]) has an invariant meaning and holds without the restriction on supp/. On the other 
hand, if supp / is small enough (but not necessarily near xq), y = ex.p^{t9) defines a local diffeo- 
morphism t0 1— > y G supp/, therefore after making the change of variables y = exp^{t6), see |SU3] . 
this becomes 

flO) N f(x)-^— [a (xv)^^^^^^^det^^dv 

(10) N^,J{x) - J Am{x,y) det ^^^^ dy, 

where 

(11) Am{x,y) = |af„(x,grad^,p(x,y)) |^, 

y are any local coordinates near supp/, and p{x,y) = |exp~^y|. Formula ([36|) can be also un- 
derstood invariantly by considering dxP and dyp as tensors. For arbitrary / G L?'{M) we use a 
partition of unity in TA/{°* to express Na^f{x) as a finite sum of integrals as above, for x near any 
fixed Xq. 

We get in particular that has the pseudolocal property, i.e., its Schwartz kernel is smooth 
outside the diagonal. As we will show below, similarly to the analysis in [SUSllSUl] . is a ^'DO 
of order —1. 

We always extend functions or tensors defined in M as outside M. Then N^f is well defined 
near M as well and remains unchanged if M is extended such that it is still in Mi, and / is kept 
fixed. 
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2.3. Decomposition of symmetric tensors. For more details about the decomposition below, 
we refer to [Shl| . Given a symmetric 2-tensor / = fij, we define the 1-tensor 6f called divergence 
of / by 

in any local coordinates, where are the covariant derivatives of the tensor /. Given an 1-tensor 
(a vector field or an 1-form) v, we denote by dv the 2-tensor called symmetric differential of v: 

[dv]ij = ^ {ViVj + VjVi) . 

Operators d and —6 are formally adjoint to each other in L'^{M). It is easy to see that for each 
smooth V with w = on 5M, we have I{dv){'j) = for any geodesic 7 with endpoints on dM. This 
follows from the identity 

(12) ^^(^y^^(^t)),m = {dvm),j\t)). 

If a = {am} is as in the Introduction, we get 

(13) la{dv) = 0, yveC^{M), 

and this can be extended to u G Hq{M) by continuity. 

It is known (see |Shl| and (jl5p below) that for g smooth enough, each symmetric tensor / G 
L^{M) admits unique orthogonal decomposition f = + dv into a solenoidal tensor Sf := and 
a potential tensor Vf := dv, such that both terms are in Lp'{M), is solenoidal, i.e., 6f^ = in 
M, and v S Hq{M) (i.e., f = on dM). In order to construct this decomposition, introduce the 
operator = 5d acting on vector fields. This operator is elliptic in M, and the Dirichlet problem 
satisfies the Lopatinskii condition. Denote by the Dirichlet realization of A'^ in M . Then 

(14) v = {^hr^5f, r = f-d{Ah)-'5f. 
Therefore, we have 

V = d{A}y)~U, S = Id-V, 

and for any g G C^(M), the maps 

(15) {A%y^ : H-\M) ^ H^{M), V,S:L'^{M) — > L^{M) 

are bounded and depend continuously on g, see |SU4l Lemma 1] that easily generalizes for manifolds. 
This admits the following easy generalization: for s = 0, 1, . . . , the resolvent above also continuously 
maps H^~^ into H'^'^^ n Hq, similarly, V and S are bounded in H'^, if 5 G C*', k ^ 1 (depending 
on s). Moreover those operators depend continuously on g. 

Notice that even when / is smooth and / = on dM, then does not need to vanish on dM. In 
particular, extended as to Mi, may not be solenoidal anymore. To stress on the dependence 
on the manifold, when needed, we will use the notation vm and as well. 

Operators S and V are orthogonal projectors. The problem about the s-injectivity of la then 
can be posed as follows: if laf = 0, show that = 0, in other words, show that la is injective on 
the subspace SL^ of solenoidal tensors. Note that by (fT3]) and ([6]), 

(16) Na = NaS = SNa, VNa = NaT = 0. 

Lemma 2. Let a = {om} with am G C^{TCm) be as in the Introduction. The following statements 
are equivalent: 

(a) la is s-injective on L'^(M); 

(h) Na : L?{M) L'^{M) is s-injective; 

(c) Na : L\M) L\Mi) is s-injective; 
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(d) IfT"^ is the set of geodesies issued from (suppom)™* as in and = UF^, then Ip" is 
s-injective. 

Proof. Let la be s-injective, and assume that Naf = in M for some / G L'^{M). Then 

= (iV./,/)L.(M)=Ell«-^/llW™,dM) =^ f'=^- 

This proves the imphcation (a) =^ (6). Next, (b) =^ (c) is immediate. Assume (c) and let / G Lp'{M) 
be such that I^f = 0. Then N^f = in Mi, therefore = 0. Therefore, (c) (a). Finally, 
(a) 44> {d) follows directly form the definition of Iq. □ 

Remark. Lemma [2] above, and Lemma [4|^a) in next section show that (suppam)™* in (d) can be 
replaced by suppom if F" is a complete set of geodesies. 

3. MiCROLOCAL PARAMETRIX OF Na 

Proposition 1. Let g = go €z C^{M) be a regular metric on M, and let Ti' T~i be as in 
Theorem 0. 

(a) Let a be as in TheoremlEfa). Then for any t = 1,2, . . . , there exists k > and a bounded 
linear operator 

Q : H\Mi)^SL\M), 

such that 

(17) QNaf = fli + Kf, MfeH\M), 

where K : H^{M) SH^+\M) extends to K : L'^{M) SH\M). Ift = oo, then k = oo. 

(b) Let a = Ug be as in Theorem\^b) . Then, for g in some neighborhood of go, (a) still holds 
and Q can be constructed so that K would depend continuously on g. 

Proof. A brief sketch of our proof is the following: We construct first a parametrix that recovers 
microlocally /|^^ from N^f- Next we will compose this parametrix with the operator i— > /|^ 
as in [SU31 ISU4j . Part (b) is based on a perturbation argument for the Fredholm equation (fT7|l . 
The need for such two step construction is due to the fact that in the definition of Z'', a solution 
to a certain boundary value problem is involved, therefore near dM, our construction is not just a 
parametrix of a certain elliptic ^'DO. This is the reason for losing one derivative in ([7]). For tensors 
of orders and 1, there is no such loss, see [SU3j and ([6T]) . ([62]) . 

As in [SU4j . we will work with ^'DOs with symbols of finite smoothness S> 1. All operations 
we are going to perform would require finitely many derivatives of the amplitude and finitely many 
seminorm estimates. In turn, this would be achieved if 5 G C^, k ^ 1 and the corresponding ^'DOs 
will depends continuously on g. 

Recall |SU31[SU4] that for simple metrics, is a ^'DO in M™* of order —1 with principal symbol 
that is not elliptic but N + \D\~^T' is elliptic. This is a consequence of the following. We will say 
that Na (and any other 'J'DO acting on symmetric tensors) is elliptic on solenoidal tensors, if for 
any (x,^), ^ 7^ 0, <Jp{Nay^''\x,^)fki = and ^Vij = imply / = 0. Then is elliptic on solenoidal 
tensors, as shown in |SU3| . That definition is motivated by the fact that the principal symbol of 5 is 
given by fij 1— > i^fiji s-injectivity is equivalent to the statement that Nf = and 5/ = in M 
imply / = 0. Note also that the principal symbol of d is given by Vj ^ {^iVj + ^jVi)/2, and o'p{N) 
vanishes on tensors represented by the r.h.s. of the latter. We will establish similar properties of 
Na below. 

Let Na^ be as in Section 12.21 with m fixed. 
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Lemma 3. N^^^ is a classical '^DO of order —1 in M{°*. It is elliptic on solenoidal tensors at 
(xo,^°) if and only if there exists 6q G T^^Mi \ with {^^,0q) = such that ao(xo,^o) 7^ 0. The 
principal symbol o'p{Na^) vanishes on tensors of the kind fij = {£,iVj +^jVi)/2 and is non-negative 
on tensors satisfying C fij = 0. 

Proof. We established the pseudolocal property already, and formulas ([9]), ([36]) together with the 
partition of unity argument following them imply that it is enough to work with x in a small 
neighborhood of a fixed xq € M™*, and with / supported there as well. Then we work in local 
coordinates near xq. To express as a pseudo-differential operator, we proceed as in |SU31[SU4] . 
with a starting point (p6]) . Recall that for x close to y we have 

p^{x,y) = Gf^{x,y){x - y)\x - yY , 

where G^]\ G^f G^f are smooth and on the diagonal. We have 

'■J *j 

G\f{x,x) = G\f{x,x) = G\f{x,x) = g,,{x). 
Then N^^ is a pseudo-differential operator with amplitude 

e-'«-^(G«z-z) ' \al{x,g-'G^'^ 

U»j (3) 

x[G(^).]jG(^).]^.[G(^).]jG(^).]^^^d., 

where G- {x,y) = G\- {y,x). As in |SU4j . we note that Mij^i is the Fourier transform of a 
positively homogeneous distribution in the z variable, of order n — 1. Therefore, Mijki itself is 
positively homogeneous of order —1 in ^. Write 

(19) M{x,y,0 = J e~'^-'\z\-''+^m{x,y,e)dz, = z/\z\, 

where 

mi,kiix,y,e) = (^G^'^e-e) ' \alix,g-'G^^^e)\' 

X [g(2) e] . 0] . [^(2) e] ^ [(5(2) e] ^ - "^^^ ^ 



z 



Y^det5(x)' 



and pass to polar coordinates z = rO. Since m is an even function of 0, smooth w.r.t. all variables, 
we get (see also jjii Theorem 7.1.24]) 



(21) M(x,y,e)=vr/ m{x,y,9)5{e-i)de. 

J\e\=i 

This proves that M is an amplitude of order —1. 

To obtain the principal symbol, we set x = y above (see also |SU31 sec. 5] to get 



(22) ap(iV^„)(x,0 = M(x,x,0=vr / m(x, x, e)5(0 • C) d0, 

J\e\=i 
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where 

(23) m'^^\x,x,e) = |af„(x,^)|^7drt^((7,j(x)e^e^)^~^e^0^'0'=^'. 

To prove ellipticity of M{x,^) on solenoidal tensors at (xo,^°), notice that for any symmetric 
real fij, we have 

(24) m}^^\x^,xo,e)fijhi = |a«„(xo,^)|Vdet5(xo) {gij{xQ)e'e^)^'\hje'e^f > 0. 

This, ([22]), and the assumption 0^(2:0, 6*0) / imply that M*-'^'(xo, xq, i^)fijfki = yields fijO'-Q^ = 
for 9 perpendicular to , and close enough to ^o- If in addition {C^Y fij = 0, then this implies 
fijO^O^ = for 9 £ neigh(0o)) and that easily implies that it vanishes for all 6. Since / is symmetric, 
this means that / = 0. 

The last statement of the lemma follows directly from (|22|). (|23|) . (|24|) . 

Finally, we note that (|23|) . (|24|) and the proof above generalizes easily for tensors of any order. □ 

We continue with the proof of Proposition [TJ Since (b) implies (a), we will prove (b) directly. 
Notice that Ti.' and 7i satisfy the properties listed in the Introduction, right before Theorem [2l if 
g = go- On the other hand, those properties are stable under small C'' perturbation of go. We will 
work here with metrics g close enough to 50- 

By Lemma [3l since T{7i') is complete, Na defined by ([6]) is elliptic on solenoidal tensors in M. 
The rest of the proof is identical to that of [SU4t Proposition 4]. We will give a brief sketch of it. 
To use the ellipticity of A^^^^ on solenoidal tensors, we complete to an elliptic ^'DO as in |SU4] . 
Set 

(25) W = Na + \D\-'VM,, 

where l-Dl""*^ is a properly supported parametrix of (—Ag)^/"^ in neigh(Mi). The resolvent (— A|^^ ^)^ 
involved in Vmi and Sm-^ can be expressed as Ri + R2, where Ri is any parametrix near Mi, and 
i?2 : -Lcomp(^i) ^ C^iMi), R2 : H\Mi) H^+^{Mi), where I = l{k) » 1, if A; > 1. Then W is 
an elliptic ^'DO inside Mi of order —1 by Lemma [3l 

Let P be a properly supported parametrix for W of finite order, i.e., P is a classical ^DO in the 
interior of Mi of order 1 with amplitude of finite smoothness, such that 

(26) PW = ld + Ki, 
and Ki : Ll^^p{Mi) H\Mi) with / as above. Then 

Pi := Sm,P 

satisfies 

(27) PiN^=Sm,+K2, 

where K2 has the same property as Ki. To see this, it is enough to apply Smi to the left and right 
of ([26]) and to use (fTGll . 

Next step is to construct an operator that recovers /|^, given and to apply it to PiA'^q — K2. 
In order to do this, it is enough first to construct a map P2 such that if /|^^ and vm^ ai'e the 
solenoidal part and the potential, respectively, corresponding to / G L'^{M) extended as zero to 
Mi\M, thenPs 

^ ^MilaA/- This is done as in [SU3J and [SU4. Proposition 4]. We also have 
P2P1 : H'^{Mi) H^''^[dM). 
Then we showed in [SU4.j Proposition 4] that one can set 

g = (id + dPP2)Pi, 
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where R : h u is the Poisson operator for the Dirichlet problem A^u = in M, u\qm = h. 

As explained above, we work with finite asymptotic expansions that require finite number of 
derivatives on the amplitudes of our ^'DOs. On the other hand, these amplitudes depend con- 
tinuously on 5 G C^, k ^ 1. As a result, all operators above depend continuously on ^ G C*^, 
A; > 1. □ 

The first part of next lemma generalizes similar results in \S\J3\ Thm 2] , \Ch.\ ISSUj to the present 
situation. The second part shows that Irf = implies that a certain /, with the same solenoidal 
projection, is flat at dM. This / is defined by the property (|29p below. 

Lemma 4. Let g G C^{M) be a regular metric, and let T he a complete set of geodesies. Then 

(a) Ker /rn5L^(M) is finite dimensional and included in C\M) with I = l{k) oo, as k ^ oo. 

(b) If Irf = with f G L'^{M), then there exists a vector field v G C\M), with v\qm = and I 
as above, such that for f := f — dv we have 

(28) d'^fldM = 0, \a\ < I, 

and in boundary normal coordinates near any point on dM we have 

(29) fni = 0, Vi. 

Proof. Part (a) follows directly from Proposition [TJ 

Without loss of generality, we may assume that Mi is defined as Mi = {x, dist(2;, M) < e}, with 
e > small enough. By Proposition [H applied to Mi, 

(30) fli, G &{Mi), 
where / ^ 1, if A; ^ 1. 

Let X = (x' , x") be boundary normal coordinates in a neighborhood of some boundary point. We 
recall how to construct v defined in M so that (|29p holds, see |SU2j for a similar argument for the 
non-linear boundary rigidity problem, and [El ISh2l ISU31 ISU4j for the present one. The condition 
(/ — dv)in = is equivalent to 

(31) VnVi + ViVn = 2fin, t;|^n=o = 0, i = l,...,n. 

Recall that ViVj = diVj — T^jVk, and that in those coordinates, = P^^ = 0. If i = n, then ([3T]) 
reduces to V„t>„ = dnVn = fnn, Vn = for = 0; we solve this by integration over < < e <C 1; 
this gives us Vn- Next, we solve the remaining linear system of n — 1 equations for i = 1, . . . , n — 1 
that is of the form VnVi = 2 fin — ViVn, or, equivalently, 

(32) d„Vi - 2r"jz;„ = 2/i„ - diVn, Vi\xr^=o = 0, z = 1, . . . , n - 1, 

(recall that a = l,...,n— 1). Clearly, if g and / are smooth enough near DM, then so is v. If we set 
f = f^ above (they both belong to Ker/r), then by (a) we get the statement about the smoothness 
of V. Since the condition (j29p has an invariant meaning, this in fact defines a construction in some 
one-sided neighborhood of dM in M. One can cut v outside that neighborhood in a smooth way 
to define v globally in M. We also note that this can be done for tensors of any order m, see [Sh2] . 
then we have to solve consecutively m ODEs. 

Let f = f — dv, where v is as above. Then / satisfies ()29p . and let 

(33) flj^=f-dvM, 

be the solenoidal projection of / in Mi. Recall that /, according to our convention, is extended 
as zero in Mi \ M that in principle, could create jumps across dM. Clearly, ff^^ = ff^^ because 
f — f = dv in M with v as in the previous paragraph, and this is also true in Mi with /, / and 
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V extended as zero (and then f = on dMi). In (|33p . the l.h.s. is smooth in Mi by (j30p . and 
/ satisfies (j29p even outside M, where it is zero. Then one can get vmi by solving (j3ip with M 
replaced by Mi, and / there replaced by flj_^ S (7'(Mi). Therefore, one gets that vmi, and therefore 
/, is smooth enough across dM, if 5 E C^, k ^ 1, which proves ([28]) . 

One can give the following alternative proof of (j28p : Let be related to F, as in Theorem [2j 
One can easily check that Na, restricted to tensors satisfying (p9]) . is elliptic for ^„ 7^ 0. Since 
^af = near M, with / extended as outside M, as above, we get that this extension cannot 
have Conor mal singularities across dM. This implies ()28p . at least when g E C°°. The case of g of 
finite smoothness can be treated by using parametrices of finite order in the conormal singularities 
calculus. □ 

4. S-INJECTIVITY FOR ANALYTIC REGULAR METRICS 

In this section, we prove Theorem [TJ Let g be an analytic regular metrics in M, and let Mi D M 
be the manifold where g is extended analytically according to Definition [H Recall that there is an 
analytic atlas in M, and dM can be assumed to be analytic, too. In other words, in this section, 
(M, dM, g) is a real analytic manifold with boundary. 

We will show first that I^f = implies E A{M). We start with interior analytic regularity. 
Below, WFa(/) stands for the analytic wave front set of /, see [Sj ITre] . 

Proposition 2. Let {xo,S,^) E T*M \ 0, and let 70 he a fixed simple geodesic through xq normal 
to . Let Lf{'y) = for some 2-tensor f E L'^{M) and all 7 E neigh(7o). Let g he analytic in 
neigh(7o) and 6f = near xq. Then 

(34) (xo,e°)0WFA(/). 

Proof. As explained in Section 12.11 without loss of generality, we can assume that 70 does not self- 
intersect. Let [/ be a tubular neighborhood of 70 with x = {x', x") analytic semigeodesic coordinates 
in it, as in the second paragraph of Section |2.1[ We can assume that xq = 0, gij{0) = Sij, and 
x' = on 7o. In those coordinates, U is given by \x'\ < e, l~ < x"^ < with some < e <C 1, and 
we can choose e <C 1 so that {x" = |x'| < e} lie outside M. Recall that the lines x' = const, in 
U are geodesies. 

Then = ((^°)', 0) with ^0 = 0. We need to show that 

(35) (O,e°)0WFA(/). 

We choose a local chart for the geodesies close to 70. Set first Z = {x" = 0; |x'| < 7e/8}, and 
denote the x' variable on Z by z' . Then z' , 0' (with \0'\ <C 1) are local coordinates in neigh (70) 
determined by {z',9') — > 7{^',o),{e',i)- Each such geodesic is assumed to be defined on l^ < t < l^, 
the same interval on which 70 is defined. 

Let Xn{z'), = 1, 2, . . . , be a sequence of smooth cut-off functions equal to 1 for \z'\ < 3e/4, 
supported in Z, and satisfying the estimates 

(36) |5"xiv| < (CAf)l"l, \a\<N, 
see [Eil Lemma 1.1]. Set 6 = {9',1), l^'l < 1, and multiply 

If (7(z',o),e) = 

by XNiz')e^'^^''^' , where A > 0, ^' is in a complex neighborhood of {S,^)', and integrate w.r.t. z' to 
get 

(37) J j e"^^'-«'x7v(^')/., b(z'fi)M %',o)A^)il',o)/*) did/ = 0. 
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For l^'l <^ 1, {z',t) £ Z X (1^,1^) are local coordinates near 70 given by x = "y[z',o),e{'t)- 
If 6' = 0, we have x = {z', t). By a perturbation argument, for 9' fixed and small enough, (t, z') 
are analytic local coordinates, depending analytically on 9' . In particular, x = (z' + t9',t) + 0{\9'\) 
but this expansion is not enough for the analysis below. Performing a change of variables in ()37p . 
we get 

(38) J e'^''^'''^'>^'aN{x,e')fij{x)b\x,9')V{x,9')dx = 

for l^'l <C 1, VA, Vi^', where, for \9'\ <^ 1, the function {x,9') 1— > is analytic and positive for x 
in a neighborhood of 70, vanishing for x ^ U, and satisfying (136p . The vector field b is analytic on 
suppaAT, and 6(0,^') = 9, a]\f{0,9') = 1. 

To clarify the arguments that follow, note that if g is Euclidean in neigh(7o), then (|38p reduces 

to 



e'^(i'^-'^'<'y-XNMx)9'0'dx = O, 



where xn = XNix' — x'^9'). Then = , —9' ■ is perpendicular to 9 = {9' , 1). This implies that 
(39) / e'^^-^XNMx)9\C)9^iC) dx = 



for any function 0(^) defined near such that 0(^)-^ = 0. This has been noticed and used before if g 
is close to the Euclidean metric (with xn = l), see e.g., |SU2] . We will assume that 9{^) is analytic. 
A simple argument (see e.g. [Shll ISU2| ) shows that a constant symmetric tensor /,-,• is uniquely 
determined by the numbers fij9^9^ for finitely many 0's (actually, for N' = (n + l)n/2 0's); and in 
any open set on the unit sphere, there are such 0's. On the other hand, / is solenoidal near xq. To 
simplify the argument, assume for a moment that / vanishes on dM and is solenoidal everywhere. 
Then CfijiO — 0- Therefore, combining this with ()39p . we need to choose = n(n — l)/2 vectors 
^(0) perpendicular to ^, that would uniquely determine the tensor / on the plane perpendicular 
to ^. To this end, it is enough to know that this choice can be made for ^ = then it would 
be true for ^ € neigh(^''). This way, Cfij{C) = and the equations ([39]) with the so chosen 
^p('?)) P = 1) • • • ) form a system with a tensor- valued symbol elliptic near ^ = The C°° ^'DO 
calculus easily implies the statement of the lemma in the C°° category, and the complex stationary 
phase method below, or the analytic ^'DO calculus in [Trej with appropriate cut-offs in implies 
the lemma in this special case {g locally Euclidean). 

We proceed with the proof in the general case. Since we will localize eventually near xq = 0, 
where g is close to the Euclidean metric, the special case above serves as a useful guideline. On the 
other hand, we work near a "long geodesic" and the lack of points conjugate to xq = along it will 
play a decisive role in order to allow us to localize near x = 0. 

Let 9{S,) be a vector analytically depending on ^ near ^ = S,^, such that 

(40) 0(0 -^ = 0, r(0 = l, 9{e) = en. 

Here and below, ej stand for the vectors d/dxK Replace 9 = {9', 1) in (j38p by 9{(^) (the requirement 
l^'l ^ 1 is fulfilled for ^ close enough to ^°), to get 

(41) [ e'^'^^^^^^dNix,Okix)b\x,^)Vix,Odx = 0, 



where oat is analytic near 70 x {C^}, and satisfies (p6]l for ^ close enough to and all x. Next, (p, 
b are analytic on suppa^v for ^ close to In particular, 

b = 7(.',o),(e'(5),i)(0' t = t{x,9'{0), z' = z'{x,9'{0), 
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and 

mo = m, a^(o,e) = i. 

The phase function is given by 

(42) ^{x,O=z'{x,0'iC)).e. 

To verify that is a non-degenerate phase in neigh(0, ^''), i.e., that detipx^{0,^^) 7^ 0, note first 
that z' = x' when x" = 0, therefore, {dz' /dx'){0, 9{S,)) = Id. On the other hand, hnearizing near 
= 0, we easily get {dz'/dx'^){0,9{C)) = -e'{(,). Therefore, 

by (l40]l . So we get ^Px^iO,£,) = Id, which proves the non-degeneracy claim above. In particular, 
we get that x 1— > if^{x,^) is a local diffeomorphism in neigh(O) for ^ G neigh(^'^), and therefore 
injective. We need however a semiglobal version of this along 70 as in the lemma below. For this 
reason we will make the following special choice of ^(^). Without loss of generality we can assume 
that 

f = e^-\ 

Set 

(43) = (6, . . . , j" + ---+^n-2 + in ^ ^ 

V ?n-l 

If n = 2, this reduces to 0(^) = (—■^2/^1; !)• Clearly, 0(^) satisfies (j^Oj) . Moreover, we have 

(44) #(e°) = e., z. = l,...,n-2, ^(e°) = 0, ^(e°) = -e„-i, 

In particular, the differential of the map 5*"^^ 9 ^ ^ ^'(0 is invertible at ^ = = e"""^. 
Lemma 5. Let 6{S) he as in ^4^; o,^^ Lp{x,(,) be as in (4^ - Then there exists 6 > such that if 

for some x G U, \y\ < 6, |^ — ^"^1 < S, complex, then y = x. 

Proof. We will study first the case y = 0, ^ = , x' = 0. Since 975(0,^) = 0, we need to show that 
(^5((0,x"),C°) = for (Cx'^) G U (i.e., for r < x" < /+) implies x" = 0. 

To compute ip^{x,i^), we need first to know dz'{x, 9')/ 89' at 9' = 0. Differentiate 7^^, ^g, ^-^{t) = 
x' w.r.t. 9', where t = t{x,9'), z' = z'{x,9'), to get 

dz' dt 
^e.7(2',o),(e',i)(0 + dzn[z'fl),{9',i)^^) • -QQ-^ + 7(2',o),{e',i)(*)^ = 0- 

Plug 6*' = 0. Since dt/d9' = at 6*' = 0, we get 



dz' 

-r^ = -5e.7(^',o),(eM)(^") 



?'=0,x'=0 



89, 

where the prime denotes the first n — 1 components, as usual; Ji^{x") is the Jacobi field along 
the geodesic x" i— > 7o(x") with initial conditions Ju{0) = 0, DJ,{0) = e,; and D stands for the 
covariant derivative along 79. Since z'{{0,x^),9'{S,^)) = 0, by (j42p we then get 



0^1 oil 
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By (jMD, (recall that = e""!), 

f-JrH^"), / = l,...,n-2, 

(45) ^((o,a;-),eO)= <^0, / = n - 1, 

where J"~^ is the (n — l)-th component of Jy. Now, assuming that the l.h.s. of (j45p vanishes 
for some fixed x" = t^, we get that J^~^(to) = 0, = 1, ... ,n — 1. On the other hand, Jy are 
orthogonal to because the initial conditions Ji/(0) = 0, DJy{0) = are orthogonal to e^, too. 
Since = ^in, this means that J" = 0. Therefore, Juito), v = l,...,n — 1, form a linearly 
dependent system of vectors, thus some non-trivial linear combination Jy{tQ) vanishes. Then the 
solution Jo(t) of the Jacobi equation along 70 with initial conditions Jo(0) = 0, £)Jo(0) = a'^Cu 
satisfies J(to) = 0. Since DJo(O) 7^ 0, Jq is not identically zero. Therefore, we get that xq = and 
X = (0, io) are conjugate points. Since 70 is a simple geodesic xq, we must have to = = 2;". 
The same proof applies if x' 7^ by shifting the x' coordinates. 

Let now y, ^ and x be as in the Lemma. The lemma is clearly true for x in the ball 5(0, ei) = 
{|x| < El}, where £1 <C 1, because 93(0, is non-degenerate. On the other hand, <f^{x, / V?(2/) C) 
forxG f7\S(0,ei), 2/ = 0, e = C°. Hence, we still have (/3g(x,^) / ^^{y,0 ^ small perturbation 
of y and ^. □ 

The arguments that follow are close to those in |KSUl Section 6]. We will apply the complex 
stationary phase method [SjJ. For x, y as in LemmaO and |?? — C^'l < (5/(7, C ^> 2, 5 <C 1, multiply 
dm) by 

where x is the characteristic function of the ball -6(0, 6) C C", and integrate w.r.t. ^ to get 

(46) / / e^"*(^'-''''«)s^(x, v)hixW{x, e)^(x, ^ dx d^ = 0. 



Here cat = ~ 'n)(^N is another amplitude, analytic and elliptic for x close to 0, 1^^ — ??| < 
and 

$ = -</'(y,0 + ^(a;,e) + ^(e - '?)^- 

We study the critical points of 1-^ If y = x, there is a unique (real) critical point .^c = "^j and 
it satisfies Q'^^^ > at = ^c- For y 7^ x, there is no real critical point by Lemma [5l On the 
other hand, again by Lemma [5l there is no (complex) critical point if |x — y| > 5/Ci with some 
Ci > 0, and there is a unique complex critical point if 1^; — y| < ^/C2, with some C2 > Ci, still 
non-degenerate if (5 <C 1. For any Co > 0, if we integrate in (|16]) for |x — y| > S/Cq, and use the 
fact that |$^| has a positive lower bound (for ^ real), we get 



(47) 



e^my,-'^,0~a^^^^ v)fij{x)b\x, 0^'(x, dx d^ 

\x-y\>5/Co 



< C3{C3N/X)^ + CNe-^/^. 



Estimate (j47j) is obtained by integrating N times by parts, using the identity 

as well as using the estimate (|36p . and the fact that on the boundary of integration in ^, the e'^* is 
exponentially small. Choose Co ^ C2. Note that 9$ > for G 9(suppx(- — ??)), and rj as above, 
as long as C 1, and by choosing Co ^ 1, we can make sure that is as close to tj, as we want. 
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To estimate ([I6|) for \x — y\ < S/Cq, set 

i{;{x,y,r]) := 

Note that = — i(y — x) + t] + 0{6), and ■ip{x,y,rj) = rj ■ [x — y) + ^\x — + 0{6). We will not 
use this to study the properties of ip, however. Instead, observe that at y = x we have 

(48) ipy{x, X, rj) = -(px{x, Tj), tpx{x, X, if) = ip^^x, if), ip{x, x, rj) = 0. 
We also get that 

(49) Q'ip{y,x,r]) >\x -y\'^/C. 

The latter can be obtained by setting h = y — x and expanding in powers of h. The stationary 
complex phase method [Sj], see Theorem 2.8 there and the remark after it, gives 

(50) [ e'^'f''^'''"^fij{x)B'^{x,a;X)dx = 0{X''/\C3N/X)^ + Ne-^/^), ViV, 

J\x-y\<S/Co 

where a = (y, ??), and i? is a classical analytic symbol [Sj] with principal part equal to b0b, up to an 
elliptic factor. The l.h.s. above is independent of A^, and choosing so that < A/(C3e) < + 1 
to conclude that the r.h.s. above is 0{e~^^^). 

In preparation for applying the characterization of an analytic wave front set through a general- 
ized FBI transform jSj], define the transform 

where, following fSj], a = (a^;, a^). It is a diffeomorphism from neigh(0, to its image, and denote 
the inverse one by a{(3). Note that this map and its inverse preserve the first (n-dimensional) 
component and change only the second one. This is equivalent to setting a = {y,ri), P = (y,C)) 
where C = ^yiv-, v)- Note that C = ^ + 0(5), and at y = 0, we have Q = rj. 
Plug a = a{f3) in ^ to get 

(51) j e'^^^'''^^fij{x)B'^{x,f3;X)dx = ©(e"^/^), 

where ijj, B are (different) functions having the same properties as above. Then 

(52) i;y{x, X, = -C> A{x, X, C) = C, ipix, x, C) = 0. 
The symbols in ()5ip satisfy 

(53) o-p(S)(0, 0, C) = 0{C) 61(C) up to an elliptic factor, 

and in particular, ap{B){0,0, (,^) = e„ e„, where dp stands for the principal symbol. 

Let 6*1 = e„, ^2, • • • > ^Af be = n(n — l)/2 unit vectors at xq = 0, normal to = e"~^ such that 
any constant symmetric 2-tensor / such that /"~^ = 0, V« (i.e., //^^ = 0) is uniquely determined 
by fij9^9^, 6 = 6p, p = 1, . . . , N . Existence of such vectors is easy to establish, as mentioned above, 
and one can also see that such a set exists in any open set in (^'')"'". We can therefore assume that 9p 
belong to a small enough neighborhood of 9i = e„ such that the geodesies [—1^, /+] 9 t i— > 7o,6»p(^) 
through xo = are all simple. Then we can rotate a bit the coordinate system such that = e"~^ 
again, and 9p = Cn, and repeat the construction above. This gives us A^ phase functions ip(p), and 
as many symbols S(p) in (fSTj) such that (1521) holds for all of them, i.e., in the coordinate system 
related to 6i = e^, we have 

(54) J e^"'^w(-./5);^^.(^)5g^(^,;3.A)dx = 0(e-"/^), p = l,...,N, 
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and by ((53]), 

(55) ap(S(p))(0,0,e°) = 0p®0p, p=l,...,N. 

Recall that 5f = near xq = 0. Let xo = Xoi^) be a smooth cutoff close enough to x = 0, equal 
to 1 in neigh(O). Integrate ^exp (iA^/^(i)(a;,/3))xo'5/ = w.r.t. x, and by (09]), after an integration 
by parts, we get 

(56) / e'^^w(-'''^xoix)h,ix)CHx,P;X)dx = 0{e-^/^), i = l,...,n, 



for Px = y small enough, where ap{C^){0,0,^^) = {£,^y. 

Now, the system of + n = (n + l)n/2 equations ([56|) can be viewed as a tensor-valued 

operator applied to the tensor /. Its symbol, an elliptic factor at (0,0,^*^), has "rows" given by 
^p^P; P = • • • ) and 6l,{S,^y , k = 1, . . . ,n. It is easy to see that it is elliptic; indeed, the latter is 
equivalent to the statement that if for some (constant) symmetric 2-tensor /, in Euclidean geometry 
(because gij{0) = 6ij), we have fijOpOp = 0, p = 1, . . . , N; and /"""^ = 0, i = 1, . . . , n, then / = 0. 
This however follows from the way we chose 9p. Therefore, (j35|) is a consequence of ([56j) . see 
[5j| Definition 6.1]. Note that in [Sj], it is required that / must be replaced by / in If 
/ is complex- valued, we could use the fact that /(3ft/) (7) = 0, and 1(9/) (7) = for 7 near 70 and 
then work with real- valued /'s only. 

Since the phase functions in (I54p depend on p, we need to explain why the characterization 
of the analytic wave front sets in [Sj] can be generalized to this vector-valued case. The needed 
modifications are as follows. We define h^^_^^{x, 13; X) = B^('py P = l,--- and /i*^-jy^^^(x, /3; A) = 

C^SI, k = 1, . . . ,n. Then {h^^^}, p = 1, . . . , N + n, is an elliptic symbol near (0, 0, ^°). In the proof 
of [Sjl Prop. 6.2], under the conditions ([^91) . ([52]) . the operator Q given by 



[QfUx,X) = // e'^('/'(^)(-'^)-'/'w(^-'3))/,^.(y,A)/i;^.(x,/3;A)dyd/3 



is a ^DO in the complex domain with an elliptic matrix-valued symbol, where we view / and Qf 
as vectors in R^"*"". Therefore, it admits a parametrix in H^^^o with a suitable tp (see fSj]). Hence, 
one can find an analytic classical matrix-valued symbol r{x,f3, A) defined near (0,0,^°), such that 
for any constant symmetric / we have 



Q(r(.,/3,A)ei"'^a)/)] = e^^^^W/, Vp 



The rest of the proof is identical to that of |Sj, Prop. 6.2] and allows us to show that (j5ip is 
preserved with a different choice of the phase functions satisfying (jl9]) . (j52|) . and elliptic amplitudes; 
in particular. 



/ 



iAV^(i)(-,/3)^2(x)/.,(x)dx = 0(e-"/^), Vi, j 



for (3 G neigh(0, ^'^) and for some standard cut-off X2 near x = 0. This proves (135]) . see (SjJ 
Definition 6.1]. 

This concludes the proof of Proposition [2l Notice that the proof works in the sane way, if / is a 
distribution valued tensor field, supported in M. □ 

Lemma 6. Under the assumptions of TheoremUl let f be such that I^f = 0. Then S A{M). 

Proof. Proposition [21 combined with the completeness of T, imply that f^ is analytic in the interior 
of M. To prove analyticity up to the boundary, we do the following. 
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We can assume that Mi \ M is defined by —£i < < 0, where x"" is a boundary normal 
coordinate. Define the manifold M1/2 D M by x" > — ei/2, more precisely, M1/2 = MU {— ei/2 < 
< 0} C Ml. 

We will show first that /l/^^^ ^ •^(^1/2)- Let us first notice, that in M1/2 /mi/2 ~ ~^''^a^i/2' 
where fMi/2 satisfies A'*t>jv/i/2 = in M1/2 \ M, v\qmi/2 ~ Therefore, vmi/2 is analytic up to 
dMi/2 in Ml/2 \ see |MN1 ISU4j . Therefore, we only need to show that Im^^^ is analytic in 
some neighborhood of M. This however follows from Proposition [2l applied to Mi/2- Note that if 
ei <C 1, simple geodesies through some x G M would have endpoints outside M1/2 as well, and by 
a compactness argument, we need finitely many such geodesies to show that Proposition [2] implies 
that i® analytic in, say, M1/4, where the latter is defined similarly to M1/2 by x" > — ei/4. 

To compare flj^^^ and = fl^, see also jSUSj |SU4], write flj^^^ = f - dvMi/2 i^i ^1/2 > and 
fii ~ f ~ d'^M in M. Then dvMi/2 ~ ~/mi/2 -^1/2 \ and is therefore analytic there, up to 
dM. Given x G 5M, integrate {dvM-i/2-''^'^) along geodesies in M1/2 \ -M, close to ones normal to 
the boundary, with initial point x and endpoints on dMi/2- Then we get that VMi/2\dM £ A{dM). 
Note that VM-1/2 ^ near dM, and taking the trace on dM is well defined, and moreover, if x" is 
a boundary normal coordinate, then neigh(O) 9 x" 1-^ VMi/ii'^^"') is continuous. Now, 

(57) flj = f- dvM = fiiy^ + in M, where w = vmi/2 - vm- 
The vector field w solves 

A^lD = 0, w\oM = VMi/2\dM ^ -^(.dM). 

Therefore, w £ ^(M), and by §7^, G ^(M). 

This completes the proof of Lemma [H □ 

Proof of TheoremUl Let Irf = 0. We can assume first that f = f^, and then / G ^(M) by 
Lemma m By Lemma [H there exists h G S~^Sf such that d^h = on dM for all a. The tensor 
field h satisfies (|29l) . i.e., = 0, Vz, in boundary normal coordinates, which is achieved by setting 
h = f — dvQ, where vq solves (j3ip near dM. Then vq, and therefore, h is analytic for small x" > 0, 
up to x" = 0. Lemma m then implies that /i = in neigh((9M). So we get that 

(58) / = dvQ, < x" < eo, with vo\x"=o = 0, 

where x" is a global normal coordinate, and < eo ^ 1- Note that the solution vq to (I58p (if exists, 
and in this case we know it does) is unique, as can be easily seen by integrating (/, 7^) along paths 
close to normal ones to dM and using (I12p . 

We show next that vq admits an analytic continuation from a neighborhood of any xi G dM 
along any path in M. 

Fix X G M. Let c(t), < t < 1 be a path in M such that c(0) = xq G dM and c(l) = x. Given 
e > 0, one can find a polygon xqXi . . . x^x consisting of geodesic segments of length not exceeding 
e, that is close enough and therefore homotopic to c. One can also assume that the first one is 
transversal to dM, and if x G dM, the last one is transversal to dM as well; and all other points of 
the polygon are in M™*. We choose e <C 1 so that there are no conjugate points on each geodesic 
segment above. We also assume that £ < Sq- Then f = dv near xqXi with v = vq hy ([58]) . As in 
the second paragraph of Section [2711 one can choose semigeodesic coordinates (x',x") near X1X2, 
and a small enough hypersurface Hi through xi given locally by x" = 0. As in Lemma HI one can 
find an analytic 1-form vi defined near X1X2, so that (/ — dui)m = 0, vi\x^=o = vo{x',0). Close 
enough to xi, we have vi = vq because vq is also a solution, and the solution is unique, see also 
(j32p . Since vi is analytic, we get that it is an analytic extension of vq along X1X2. Since / and 
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vi are both analytic in neigh(xiX2), and / = dvi near xi, this is also true in neigh(xiX2). So we 
extended vq along xoXiX2, let us call this extension v. Then we do the same thing near X2X3, etc., 
until we reach neigh(x), and then f = dv there. 

This defines v in neigh(x), where x € M was chosen arbitrary. It remains to show that this 
definition is independent of the choice of the path. Choose another path that connects some 
yi £ dM and x. Combine them both to get a path that connects xi € dM and yi G dM. It 
suffices to prove that the analytic continuation of vq from xi to yi equals vq again. Let ci U 71 U 
C2 U 72 U • • • U 7fc U Cfc+i be the polygon homotopic to the path above. Analytic continuation along ci 
coincides with vq again by ([58l) . Next, let pi, p2 be the initial and the endpoint of 71, respectively, 
where pi is also the endpoint of ci. We continue analytically vq from neigh(pi) to neigh(|?2) along 
7i, let us call this continuation v. By what we showed above, f = dv near 71. Since //(71) = 0, 
and v{pi) = 0, we get by (fT2]) . that {v{p2),'yi{l)) = as well, where / is such 7i(/) = P2- Using 
the assumption that 71 is transversal to dM at both ends, one can perturb the tangent vector 
7i(/) and this will define a new geodesic through p2 that hits dM transversely again near pi, where 
V = Vq = 0. Since T is open, integral of / over this geodesic vanishes again, therefore (^(^2)5 ^2) = 
for ^2 in an open set. Hence v{p2) = 0. Choose q2 € dM close enough to p2, and r]2 close enough 
to .^2 (in a fixed chart). Then the geodesic through (925 will hit dM transversally close to pi, 
and we can repeat the same arguments. We therefore showed that v = on dM near p2. On the 
other hand, vq has the same property. Since f = dv = dvo there, by the remark after ([581) . we get 
that V = Vq near p2. We repeat this along all the legs of the polygon until we get that the analytic 
continuation v of vq along the polygon, from xi to yi, equals vq again. 

As a consequence of this, we get that f = dv in M with v = on dM. Since / = f^, this implies 

/ = o. 

This completes the proof of Theorem [TJ □ 



5. Proof of Theorems [2] and [3] 

Proof of TheoremlM Theorem [2|^b), that also implies (a), is a consequence of Proposition [H as 
shown in |SU4| . see the proof of Theorem 2 and Proposition 4 there. Part (a) only follows more 
directly from [Tall Prop. V.3.1] and its generalization, see |SU31 Thm 2]. □ 

Proof of Theorem\^ First, note that for any analytic metric in Ir^ is s-injective by Theorem [TJ 
We build Qs as a small enough neighborhood of the analytic metrics in Q. Then Qg is dense in Q 
(in the C''{Mi) topology) since it includes the analytic metrics. To complete the definition of Gs, 
fix an analytic gQ £ G- By Lemma [U one can find TC' (sTi related to g = gQ and Tg, satisfying the 
assumptions of Theorem [21 and they have the properties required for g close enough to gQ. 

Let a be as in Theorem [2] with a = 1 on TC' . Then, by Theorem [2], Ia,g is s-injective for g 
close enough to gQ in C'^(Mi). By Lemma[2l for any such g, 1^^ is s-injective, where = F('H"), 
= suppo. If g is close enough to C F^ because when g = go, F" C T(Ti.) d Tg^, 

and Tg depends continuously on g in the sense described before the formulation of Theorem [3l 
Those arguments show that there is a neighborhood of each analytic go £ G with an s-injective /pg- 
Therefore, one can choose an open dense subset Gs of G with the same property. □ 

Proof of Corollary [TJ It is enough to notice that the set of all simple geodesies related to g depends 
continuously on g in the sense of Theorem [3j Then the proof follows from the paragraph above. □ 
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6. X-RAY TRANSFORM OF FUNCTIONS AND 1-FORMS/ VECTOR FIELDS 

If / is a vector field on M, that we identify with an 1-form, then its X-ray transform is defined 
quite similarly to ([T]) by 

(59) /r/(7)= r(/(7W),7(t))dt, 7 e T. 

Jo 

If / is a function on M, then we set 

(60) /r/(7)= t' f{i{t))dt, 7er. 

Jo 

The latter case is a partial case of the X-ray transform of 2-tensors; indeed, if / = ag, where / is a 
2-tensor, a is a function, and g is the metric, then /p/ = -^po, where in the l.h.s., Ip is as in ([1]), and 
on the right. It is as in (160p . The proofs for the X-ray transform of functions are simpler, however, 
and in particular, there is no loss of derivatives in the estimate ([7]), as in |SU3j . This is also true 
for the X-ray transform of vector fields and the proofs are more transparent than those for tensors 
of order 2 (or higher). Without going into details (see |SU3j for the case of simple manifolds), we 
note that the main theorems in the Introduction remain true. In case of 1-forms, estimate ([7]) can 
be improved to 

(61) \\n\LHM)/C < WN^fWmiM,) < C||/1Il2(m), 
while in case of functions, we have 

(62) \\f\\LHM)/C < WN^fWmiM,) < C||/||l2(m). 

If (M, dM) is simple, then the full X-ray transform of functions and 1-forms (over all geodesies) is 
injective, respectively s-injective, see |Mu2l i MuRl IBGl I ARj . 
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